Ballistic interferences in suspended graphene 
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Graphene is a 2-dimensional (2D) carbon allotrope with the atoms arranged 
in a honeycomb lattice |jT|. The low-energy electronic excitations in this 2D 
crystal are described by massless Dirac fermions that have a linear dispersion 
relation similar to photons [2^ ^ . Taking advantage of this "optics-like" electron 
dynamics, generic optical elements like lenses, beam splitters and wave guides 
have been proposed for electrons in engineered ballistic graphene |4l,i5j. Tuning 
of these elements relies on the ability to adjust the carrier concentration in de- 
fined areas, including the possibility to create bipolar regions of opposite charge 
{p-n regions). However, the combination of ballistic transport and complex 
electrostatic gating remains challenging. Here, we report on the fabrication and 
characterization of fully suspended graphene p-n junctions. By local electro- 
static gating, resonant cavities can be defined, leading to complex Fabry-Perot 
interference patterns in the unipolar and the bipolar regime. The amplitude 
of the observed conductance oscillations accounts for quantum interference of 
electrons that propagate ballistically over long distances exceeding 1 yum. We 
also demonstrate that the visibility of the interference pattern is enhanced by 
Klein collimation at the p-n interface |[6|, [7] . This finding paves the way to more 
complex gate-controlled ballistic graphene devices and brings electron optics in 
graphene closer to reality. 
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FIG. 1. a, Low-energy dispersion in pristine graphene. b, Sketch of an ideahzed graphene Fabry- 
Perot resonator, c, Normahzed conductance through an ideal graphene cavity as a function of the 
dimensionless parameter kpL/ir for three different finesses: F = 1600 in yellow, F = 4.93 in purple 
and F = 0.39 in blue. 



The wave nature of any physical entity, be it light or electrons, is demonstrated by 



measuring an interference pattern in an apparatus, called interferometer. We distinguish 
two types of interferometers: two path and multi-pass interferometers. Examples for the 
former are the Michelson and the Mach-Zehnder interferometer, which are both widely used 
in optical experiments. In contrast, the Fabry-Perot (FP) interferometer is of the multi-pass 
type. In optics an FP interferometer is an element that consists of a transparent glass plate 
with two partially reflective surfaces on opposite sides. It can therefore be viewed as a cavity. 
Light within the cavity is bouncing back and forth between the mirrors, but at each reflection 
a fraction of light is coupled out. The superposition of all outgoing waves gives rise to an 
intensity modulation (the interference pattern) that depends on the wavelength of the light 
and the distance between the mirrors. In solid-state physics a two-dimensional conducting 
sheet between two electrodes could realize an electronic FP interferometer. However, the 
mean free path (mfp) and phase-coherence length are usually much shorter than the device 
size leading to purely diffusive incoherent electronic transport. Before the advent of ID 
and 2D carbon allotropes, long mfps and optic-like experiments have only been realized in 
complex semiconductor heterostructures E] • During the last decade the long mfp (several 
hundreds of nm) accessible in low- dimensional carbon materials enabled one to explore 
different ballistic phenomena. In particular, FP interferences have been demonstrated in 
carbon nanotubes [10] and graphene nanostructures [7j. In the former case, interferences 
were observed over a distance as long as 0.5 /xm, whereas the cavity size was less than 100 nm 
in the latter. Nowadays, graphene junctions show a mfp close to l/im [TT] and electron 
guiding [12] as well as transverse magnetic focusing have been revealed [13] . However, even 
if graphene appears as a perfect platform for electron optics [U [5] , quantum interferences 
with long mfp > l/im and complex electrostatic gating still remain to be demonstrated. 

In this letter we report FP interferences in a gate-tunable graphene p-n junction with 
coherent ballistic electron motion. We observe an interference pattern in both the unipolar 
and bipolar regime and for all combinations of cavity mirrors, either defined at the source 
and drain contacts or at a p-n junction in the middle of the device. 

In graphene, charge carriers with small quasi-momentum Kk counted from the corner 
of the Brillouin zone have a linear dispersion. In the conduction and valence band the 
energy -E'(k) as a function of wavevector k is given by -Ec(k) = hvpk and -Et,(k) = —hvpk, 
respectively, where fc = |k| and vp is the Fermi velocity; see Fig. [Ik. These dispersion 
relations yield an electron number density n given by n = sgn[E)k'^/TT. Taking a ballistic 



graphene cavity, as sketched in Fig.[T]D, the phase difference between two successive outgoing 
rays is A0 = 2kLcos{d), where L is the size of the cavity and 9 the angle of incidence of 
the electron waves relative to the normal of the mirror. Summing over all partial outgoings 
waves, one obtains for the transmission probability |14j : 

^^^^ = l + Fsin^(A0/2) ' ^^^ 

where F = 4|ri||r2|/|tip|t2p with |ri^2| and |ti,2| the reflection and transmission amplitudes 
at the two interfaces. F is a measure of the quality factor and known as the finesse. The 
finesse of an FP interferometer determines the visibility of the interference pattern, i.e., 
the difference in T between constructive and destructive interference relative to the mean 
transmission probability. The visibility is very close to unity if F ^ 1 and it is small and 
given by F in the opposite limit. Since the electric measurement of the graphene conductance 
cannot distinguish the direction of the waves within the cavity, one has to integrate over the 
angles to deduce the expected visibility for the conductance (Gmax — G'min)/(G'max + Cmin)- 
The normalized dimensionless conductance Gn = (V^) /^(^) cos{9)d9 is shown in Fig. uJp 
as a function of the dimensionless parameter kL/n for three cases (ri = r2 = 0.8, F = 1600), 
(ri = r2 = 0.4, F = 4.9) and (ri = r2 = 0.2, F = 0.4). For all three cases resonances appear 
equidistantly whenever kpL is a multiple of vr, indicating that the conductance peaks are 
mainly due to electrons at small incident angles (see Supplementary Material for details). 

We now focus on a graphene FP cavity that consists of two segments: left and right. 
Assume that in both segments the carrier density and its sign (n and p) can be controlled 
independently by two gate voltages V^eft and Kight with identical efficiency. Figure [2] shows 
in a schematic 2D conductance map the expected FP interference patterns as a function of 
Vieft and Kight- This map has four quadrants corresponding to the four polarity configura- 
tions [15]. If the left and right segments have the same polarity {nn or pp), the cavity is in 
the unipolar regime. In contrast, in the case of opposite polarities {np or pn), the cavity 
is in the bipolar regime. In the unipolar regime only interferences of electrons bouncing 
back and forth between the two outer contacts are expected (lower left and upper right 
panel in Fig. |2]). The corresponding beating pattern evolves along Vieft = Kight and is indi- 
cated by the light-blue lines in the central panel of Fig. |2j In the bipolar case, assuming a 
semi-transparent p-n interface between the segments, the junction can be considered as two 
cavities in series, each controlled by its respective gate. Resonances in the left (right) cavity 
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FIG. 2. Schematic representation of all possible resonance conditions (central panel) in a FP 
resonator that consists of two segments of equal size (left and right panels). The carrier density 
and carrier type (n or p) in both segments can be tuned independently by V^eft and Kight- FP 
resonances may occur due to reflection at the outer contacts or at an internal n-p (p-n) interface 
between the two segments. The possible resonant conditions are indicated by color-coded lines for 
the respective cases in the middle part. 



depend only on Vmi (K-ight) and are represented by green (blue) lines. Another distinct set 
of resonances in the bipolar regime can arise from charge carriers that tunnel through the 
central p-n interface in the middle and bounce back and forth between the outer contacts. 
The corresponding beating pattern should evolve along the condition fcieft ~ bright- In terms 
of gate voltages this condition corresponds to VJeft ~ — V^ight as indicated in Fig. [2] by the 
orange lines. 

To realize experimentally such a tunable FP cavity, devices were prepared by combining 
a mechanical transfer process [16] with a hydrofluoric acid-free suspension method proposed 
by Tombros et al. [17J. The fabrication process is schematically shown in Fig. Isk. For 
a complete description, see Methods. Clean graphene is obtained by an in-situ current 
annealing process. Figure |3)d shows a colored optical image of several junctions realized 
within the same graphene flake indicated by white dashed lines. An SEM micrograph of 
three p-n junctions is presented in Fig. ^. A voltage applied on the bottom gate (Vright, 
yellow electrodes in Fig. |3}3,c) will tune the charge carrier density on the right side of a 
junction, while a voltage applied on the back gate (Vieft) will act on the left side. The device 
was measured in a ^He cryostat at a temperature of T = 1.5 K. Differential conductance, 
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FIG. 3. a, 3-diniensional illustration of the device fabrication consisting of the following steps: 
1) bottom gates are realized by standard e-beam lithography on a doped Si wafer with 300 nm 
thick Si02 top layer and 2) covered by a 700 nm thick LOR (lift-off resist) layer. 3) Graphene is 
exfoliated on a separated wafer covered by a stack of PVA/PMMA, then aligned to the bottom gates 
and transferred. 4) Contacts to the graphene are then realized with standard e-beam lithography 
followed by an evaporation of 50 nm palladium. 5) Finally, the device is suspended by selectively 
exposing the LOR with a large dose and developing it. b, Optical image of several junctions defined 
within the same graphene flake indicated by white dashed lines. Bottom gates are in yellow and 
ohmic contacts in blue. The LOR resist appears in dark green. Scale bar, 5 /im. c. Scanning- 
electron microscopy image of three graphene p-n junctions in a row. Bottom gates are shown in 
yellow and contacts in light blue. Scale bar: 1.2 //m. 



G = dl/dV was measured by standard lock-in techniques with a voltage bias of 100 /iV at 
77 Hz. Series and contact resistances were not subtracted from the data. 



Figure H^ shows the main experimental result of this work. In part a) the differential 
conductance between the two contacts as a function of Kight and V^eft is presented for a 
junction of length L = 1.2 /im and width W = 3.2 /im. As expected, this 2D plot reveals 
four regions indicated by the labels pp, nn, np and pn corresponding to different carrier 
types in the two sides of the sample. The borders between the unipolar and bipolar regions 
coincide with the charge neutrality point of the left and right graphene segments (indicated 
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FIG. 4. a, Two terminal differential conductance G = dl/dV as a function of Vright and Vjeft 
at T = 1.5 K, demonstrating independent control of carrier type and density in the left and the 
right side of the graphene sheet. Labels in each of the four conduction regimes indicate the carrier 
type, b, Conductance as a function of charge carrier density n along the dashed white line of a. 
c, Zoom into the transconductance map presented, in d, around the charge neutrality point, d, 
Transconductance obtained from the numerical derivative with respect to Kight of a emphasizing 
conductance oscillations. Color-coded lines link the observed conductance oscillations with the FP 
interference pattern expected from the simple picture of Fig. [2j light blue lines in the unipolar 
regime, green and blue lines for interferences on either side of the graphene in the bipolar case and 
orange lines for full length interference in the bipolar regime. 



by blue and green arrows respectively). The border lines between the different regions are 
not perpendicular to each other, as expected for independent gates and discussed above in 
Fig. |2]for two reasons. First, both the back gate and the bottom gate are at a distance to 
the graphene sheet that is comparable to the length of the device. This leads to a cross 
coupling between the two gates i.e., the right gate also affects the left graphene segment 
and vice versa. Second, the back gate (left gate) is below the bottom gates and therefore 
strongly screened by these, leading to a much weaker gate coupling from the back gate to 
the graphene than from the bottom gate. 

A slice through the 2D conductance plot along the white dashed line is shown in Fig. |4)d 
as a function of carrier density n which we estimated from a parallel plate capacitor model. 
A positive sign refers to the n-region and a negative sign likewise to the ]?-region. A change 



in gate-voltage of IV along the white dashed line induces a change in density of n ?^ 9 x 
10^ cm~^. The graph exhibits a remarkably sharp conductance dip, corresponding to the 
charge neutrality in the entire graphene sheet. The fact that this dip also occurs close to zero 
gate-voltage reflects the high quality of the graphene sample. The conductance saturation 
and the asymmetry at high carrier density suggest an n-type doping of the graphene contacts 
below the palladium contacts [IHIIIS]. By defining the mobility as /x = 1/exda/dn with a the 
graphene conductivity, we estimated a mobility of yU ~ 150 x 10^ cm^/Vs at ra = 4 x 10^ cm~^. 
The junction reaches a minimal conductivity of cTmin ~ 2e^ /h close to the ballistic limit of 
Ae'^/irh [20] indicating a weakly disordered graphene sheet [TT| m\. 

The most striking features in the conductance map are, however, the oscillation patterns 
visible in both the unipolar and bipolar configuration. To emphasize these features, the 
transconductance dG/dVrigu numerically calculated from the data in Fig. Ilk is presented 
in Fig. |4ji. The different interference patterns expected from the simple picture of Fig. [2] 
are well revealed in this map. In the unipolar case, one global interference pattern is vis- 
ible (highlighted in a small region by light-blue lines) corresponding to FP interferences 
of electron waves that bounce back and forth between the outer metallic contacts. In the 
data the density separation An ^ 10^° cm~^ between constructive interference peaks at 
around n = 10^^ cm~^. Based on the equation n = sgn{E)k'^/TT and the resonance condition 
AkLcavity = TT we dcrivc the relation An = 2-\/7m/Lcavity This yields an effective cavity 
length Lcavity ~ 1.1 /im in good agreement with the geometrical length L = 1.2 /im of the 
junction . In the bipolar regime at high doping two clear oscillation patterns are visible. One 
pattern is mainly tuned by K-ight, indicating that resonances arise from the cavity created 
at the right side of the p-n interface. This pattern is highlighted by blue lines in Fig. |4]i. 
The other pattern marked by green lines is predominantly tuned by V^eft, indicating that 
it originates from the cavity formed at the left side of the p-n interface. The area in the 
transconductance map marked by the dotted square is depicted enlarged in Fig. ^. It shows 
yet another oscillation pattern (orange lines). It follows gates values such that nieft ~ — "^iright 
with an interference period An ^ 3 ■ 10^ cm^^ a.t n = 10^° cm^^ leading to Lcavity ~ 1-2 fj^m. 
We conclude that these oscillations must arise from quantum interference over the full device 
length, i.e. charge carriers that tunnel through the p-n interface, changing their character 
from n to p and back. 

We emphasize at this point that there is no evidence of disorder in the transconductance 
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plot of Fig. |4]i. All patterns are regular and can be ascribed to distinct FP cavities, some of 
which extend in length over the whole device. The occurrence of FP interference patterns 
implies that the phase-coherence length exceeds twice the system size. Random disorder 
would then generate further interference patterns, but these are never regular. The fact 
that we can deduce from the interference pattern effective cavity-lengths from contact to 
contact and that we do not see any irregular features are strong signs for ballistic coherent 
transport over distances exceeding 2 fun. 
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FIG. 5. a, Schematic view of charge carrier configuration at the unipolar point (Vlcft, Kight) = 
(— 15,— 10)V indicated by a blue square in d. b. Calculated charge carrier profile at this point, 
and c, associated simulated transmission function versus incident angle 9. d, 2D conductance map 
obtained by tight-binding simulation as a function of Vieft and Kight- e,f,g. Same as a,b,c but for 
parameters at the bipolar point (Meft, ^ight) = (—15, 2)V indicated by a green square in d. 



To compare the experiment with theory, the charge transport through the graphene device 
was numerically calculated in the phasecoherent ballistic regime using the actual geomet- 
rical parameters of the device (see Supplementary Material). The computed conductance 
is shown in Fig. [5]i, which presents a purely ballistic tight-binding transport simulation 
using a real-space Green's function method within the Landauer-Biittiker formulation [22] 
combined with the quantum capacitance model for carrier density computation [23]. Since 
no dephasing or energy averaging has been incorporated, the obtained conductance modu- 



lations show the maximal possible visibility for this device geometry. The correspondence 
between theory and experiment is remarkably good and supports the picture of ballistic 
motion drawn from experiment. 

As pointed out before, we do see all the expected FP resonance patterns due to the 
formation of distinct cavities having different doping and different types of mirrors. However, 
the strength of the oscillation pattern (the visibility) varies in different regions. This is linked 
to the reflection properties of our cavity "mirrors", determining the cavity finesse. Cavities 
can be created either by sharp or smooth unipolar or bipolar potential steps. 

Let us first consider a unipolar potential step. If d is the length over which the carrier 
density varies and kp the mean Fermi wavevector, the junction is assumed smooth if kpd ^ 1 
and sharp in the opposite limit. In the smooth case the electrons follow adiabatically the 
potential evolution. They may bend similar to light rays in a medium with a slowly varying 
refractive index, but they will not backscatter. In contrast, for a sharp potential step one 
might expect backscattering. However, charge carriers in graphene are chiral, i.e. they have a 
pseudospin that is intimately tied to the wavevector. Scattering from a state with wavevector 
k into — k is forbidden, because this would flip the pseudospin. The conservation of pseudo- 
spin is only lifted, if the potential step is abrupt on the atomic scale |21] , which is not the 
case here. Thus pseudospin conservation implies that effectively all unipolar potential steps 
are highly transmissive for almost every incident angle [25] (only electrons with incident 
angles close to 7r/2 have an appreciable backscattering probability). Cavities created by 
unipolar interfaces have consequently a small finesse and will display small visibilities. 

We turn now to the reflection properties of a p-n junction. We assume that the carrier 
densities on both sides are of opposite sign but equal in magnitude and that kp denotes 
the wavevector in the homogeneous region outside the junction. In the case of a sharp 
p-n interface with kpd <^ 1 the transmission probability is given by T{6) = cos^{6) [2^ . 
The corresponding cavities have a small finesse for electron waves over a large range of 
incident angles around ^ = and will therefore have a small visibility. The situation 
is different for a smooth p-n junction. Since the charge carrier density evolves smoothly 
through zero, essentially all electron trajectories are adiabatically refracted off the interface 
(total reflection) except for the normal one, which can Klein-tunnel through the junction [6l 
[21]. The transmission probability through such a smooth p-n junction is given by the 
expression T{9) = exp(— Tr/cpcisin^^) [21] which also predicts perfect transmission at normal 
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incidence. Similar to the sharp interface one could expect the visibility of the interference 
pattern to be small. However, Klein tunneling does not only yield full transmission at 
normal incidence but also leads to a strong collimation in transmission. Therefore, almost 
all trajectories that are incident under a small angle {6 ^ 0) are exponentially suppressed and 
do not add to the total conductance. This suppression of the background current increases 
the overall visibility. Although the finesse vanishes for ^ = the contribution from small 
angles (with 6 y^ 0) yield a higher visibility than is the case for unipolar or sharp bipolar 
cavities (see Supplementary Material for further details). 

Based on the preceding discussion, we now compare the observed visibilities for two 
particular cases in a) the unipolar and b) bipolar regime at gate voltages (Vjeft, Kight) = 
(—15, — 10)V and (—15, 2)V, respectively. These two points are marked in Fig. [sli by light- 
blue and green squares. For both cases, the reflection properties at the contacts are im- 
portant parameters and need to be known. The full comparison of the simulation with the 
experiment reveals that the contacts dope graphene as n-type. In addition to this doping the 
simulation includes a small mass term in the contacted region p6]. For the unipolar case, the 
simulation yields the charge carrier density proflle plotted in Fig. [5]d. In this conflguration 
the interior of the FP cavity is p-tjpe with mirrors created by the p-n interfaces close to 
the Pd contacts. The latter can be considered to be relatively sharp due to the rapid decay 
of the contact-induced screening potential [2S]- Consequently, these interfaces are almost 
transparent and give rise to low flnesse for the FP interferences as found in the experimental 
and the numerical conductance map with respective average visibilities of 0.5% and 2%. 

In the bipolar case the simulation yields the carrier density proflle plotted in Fig. |5]f. 
It presents an extremely smooth p-n interface created electrostatically in the middle of the 
graphene sheet. The two cavities created on the right and left side will present a much higher 
flnesse than in the unipolar conflguration which consistently leads to higher visibilities in 
both the experimental and numerical conductance maps. The average visibility for the ex- 
periment is around 5%, whereas the simulation yields 20%. The calculated angle dependent 
transmission probabilities for both regimes are plotted in Figs. ^ and g respectively, which 
clearly demonstrate that the visibility difference between the unipolar and bipolar regime is 
a direct consequence of the exponential collimation of Klein tunneling [HI [Ml 125] • 

In conclusion, we have fabricated and characterized a suspended graphene p-n junc- 
tion. Fabry-Perot type conductance oscillation patterns visible in both unipolar and bipolar 
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regimes point to an extremely long mean free path and phase-coherence length > 2 /xm 
accessible in this device. Moreover we have shown that the visibility difference between 
the unipolar and the bipolar regime is due to Klein collimation occurring at a smooth p-n 
interface in graphene. Our experiment paves the way to more complex engineered ballistic 
graphene to realize electron lensing, focusing P] and efficient guiding [51 fT2]. 
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METHODS 

Suspended graphene p-n junctions were obtained by combining a mechanical transfer 
technique with a hydrofluoridric acid free suspension method developed by Tombros et 
al. [T7]. The fabrication starts with the spinning of a polymer stack consisting of a water- 
soluble PVA and PMMA layer on an oxidized (300 nm) Si wafer, on which alignment markers 
are predefined. Graphene sheets are then mechanically exfoliated on this stack and located 
by optical means. Depending on the location of the selected graphene fiake, 5/45 nm thick 
Ti/Au bottom gate electrodes are tailor-made on a separate oxidized Si wafer by standard e- 
beam lithography and covered by a 1/im thick LOR resist. Fig. [3^-1,2. Afterwards, the wafer 
that contains the graphene is immersed in deionized water. Once PVA is dissolved, the wafer 
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sinks to the bottom, leaving the PMMA-graphene membrane floating on top at the hquid-air 
interface. Next, this membrane is fished out by a transfer glass, dried and fixed in place of 
the UV-mask in a modified MJB3 mask aligner. The other wafer, containing the bottom 
gates and LOR layer, is fixed at the wafer side of the mask aligner. The graphene fiake 
is then carefully aligned to the bottom gates and finally transferred in contact by heating 
the stack to 120 °C, Fig. [3^-3. Once transferred, PMMA is dissolved in xylene. Palladium 
contacts (50 nm in thickness) are then fabricated through PMMA e-beam lithography with 
xylene as developer and lift-off solvent. Fig. [3^-4 and optical image in Fig. [3]d. In the final 
step the LOR layer below the graphene is e-beam exposed and developed in ethyl-lactate 
to obtain suspended graphene junctions. Fig. [3^-5. After mounting the finished device in a 
variable temperature ^He cryostat, we perform current annealing by applying a DC current 
at 1.5K. In this procedure the current is increased until the conductance as a function of a 
gate voltage V^ate shows a pronounced dip around Vgate = OV, reminiscent of a high quality 
device. This requires current densities that usually reach values up to 350/iA//im. 
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